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The discovery of superconductivity in A,Fe».,Se» (A = K, Cs, Rb, TI/Rb, TI/K) FS] triggered 
another surge of research on iron-based superconductors, which previously were only featured 
by iron pnictides ^?! and B-FeSe "7. The new series of superconductors can be regarded as 
formation from the intercalation of metals between the FeSe layers of B-FeSe. In comparison 
to the iron pnictides, metal-intercalated iron selenides are much complicated in terms of 
structure, chemical composition and phases. The nature of the superconducting (SC) phase, 
for example, is still in debate though considerable progress has been made over the last a few 
years [14-16], 

One of the most controversial issues is whether the SC phase is Fe-vacancy-free or not, 
i.e. the FeSe layers in A,Fe».ySey are stoichiometric (y ~ 0); or off-stoichiometric due to the 
existence of considerable Fe vacancies. The origin of this issue is largely due to the phase 
separation that inevitably occurs in these systems at 500 ~ 578 K l”, leading to the 
coexistence of the insulating phase (A2Fe4Se5) and the SC phase [1823] The thus-obtained SC 


phase is not standing freely, instead it intergrows with the insulating phase in the form of 


nanostrip and its volume fraction is quite low, 10 - 20%, as estimated by various 


measurements [°*], This is the main obstacle that prohibits the precise determination of the 
structure and composition of the SC phase. Ying and coworkers 0. 311 verified that the SC 
phases in the K-intercalated iron selenides are almost no Fe vacancy in the FeSe layers based 
on their study of the superconductors obtained via a liquid ammonia route. The neutron 
diffraction showed that the 43 K SC phase in Li/ammonia co-intercalated FeSe compound is 
Fe-vacancy-free in a separate study P”. Guo et al. P?! synthesized a 37 K SC phase containing 
slight Fe vacancy by the similar liquid ammonia method. In the K,Fe» Se» film grown on 
SrTiO; (001) substrate, an Fe-vacancy-free phase KFe.Se> with N2 x V5 charge ordering was 
observed in the SC region 134, 35] Different opinions, however, exist. For example, the SC 
phases are thought to originate from superstructures due to Fe vacancy, such as 2x 4 D6) or N8 
x V10 P7. Meanwhile, phases with disordered Fe vacancies were also reported to be SC PS 31, 
Apart from Fe vacancies, the possibility of the existence of about 2096 excess Fe in an SC 
phase with transition temperature (T.) of 44 K in the K-Fe-Se system was also reported PL 

Another puzzling issue is the existence of several SC phases with different T. from 30 K 
to 46 K, which are often mixed in a single sample. For instance, apart from a dominated 30 K 
phase, a 44 K phase of trace amount sometimes can be observed in a K,Fe»,Se» sample 
obtained by high-temperature route ll The amount of 44 K phase can be enhanced in 
samples with a little less K content (x ~ 0.6 - 0.7), but its real compositions are unknown. The 
results in ref. [31] revealed that the obtained SC phases differ only in K contents, Ko.3Fe2Se2 
with a T. = 44 K and KogsFe5Se» with a Te = 30 K. Superconductivity with T; of 30 K and 43 
K was also observed in K,Fe» Se» samples with ND x Ds] superstructure due to K vacancy 
ordering, but the phases for different T, have not been specified "^ ?* “* 9! In addition, K- 
intercalated iron selenide with excess Fe was also proposed to have T. of 44 K 10], 

So far, there has been no consensus on the answers to the above issues. Moreover, the 


underlying mechanism behind the lattice stability and vacancy are also poorly understood, in 


particular, the origin for multifarious structures. To address these issues, here, we first study 
the energetic change and structural evolution as a function of intercalated metal content by 
taking K,Fe>.,Se> as an example using first-principles calculations. Two competing factors are 
found to dominate the formation of the phases and structural evolution. One is due to the 
energy increase caused by the accumulation of negative charge in FeSe layers. The other is 
due to the Coulomb attraction between K ion layers and negatively charged FeSe layers. Then 
we show that the intercalated K content at 0.25 € x < 0.6 can stabilize the body centered 
tetragonal (bct) lattice and Fe is favored to fully occupy its site. At x > 0.6, the structure is 
stabilized by creating Fe vacancies in FeSe layers. At low intercalated K level, x « 0.25, the 
structure will collapse due to the lattice instability. Then a schematic phase diagram is 
constructed accordingly whereby we speculate the SC phases are K,Fe» Se» with 0.25 <x < 
0.6 and y « 0. Our finding sheds some light on understanding this distinct SC family and 
should be also applicable to other metal-intercalated iron selenides besides K,Fe» Se». 

To begin with, we study the formation energy for the process of K intercalation. The 
process can be expressed as a chemical reaction: 

xK + 2FeSe — K,Fe2Se2 (0 < x x 1) 

The formation energy per unit cell can be described as AE; = Exxrese? - Exore2se2 - XEx, where 
Exxre2se2 and Exore2se2 are the total energies of KxFe2Se2 and an assumed KoFe2Se2 with a 
similar bct structure but without any K ion, Ex is the energy of elemental K. The variation of 
AE, vs. K content is shown in Figure 1, which indicates that the intercalation of K into FeSe 
layers is always energetically favored for the bct structure. 

Then we break down AE; in terms of energy from the structural constituent units and 
their interactions. For K intercalation process, as shown in Figure 1(a), K atoms lose their 
valence electrons and form K ion layers between adjacent FeSe layers after K entering into 


the lattice, while FeSe layers are charged and deformed, and eventually leading to K,Fe»Se» 


with a bct structure. Hence, the following contributions to the total formation energy are 
considered: 

1) The formation of K ion layers: AEk ion layers = Eck ion layers) - XEK; 

2) The electron doping in FeSe layers: AE.. doping = EiFese layers)x- - EFeSe layers} 

3) The deformation of FeSe layers: AErese layer deformation = Efese layers - EKore2se2; 

4) The Coulomb attraction between K ion layers and FeSe layers: AEc = Ekxrose - Eirese 
layers)x- ~ Eck ion layers)x+- 

And Ek ion layers)x+ 18 the energy of K ion layers, E(Fese layers)x- the energy of charged FeSe layers, 
and Erese layers the energy of FeSe layers in K,Fe2Se2. Therefore, the total formation energy 
can be written as the sum of these four energy changes: AE; = AE ion layers + AEe- doping + 
AErese layer deformation + AEc. 

The variations of these four energy contributions as a function of K content are 
calculated and the results are shown in Figure 1(b), respectively. First of all, FeSe layers are 
prominent in contributing more and more positive energy to AE; when they are negatively 
charged with more and more electrons. This is easily understood since charging a neutral 
FeSe sheet will incur additional energy like charging a capacitor. Inversely, with the 
increasing K content, more charges will come into effect in the Coulomb attraction between K 
ion layers and the charged FeSe layers and thus greatly enhance |AEc|, which will offset 
sufficiently the increase of energy induced by electron doping into the FeSe layer. In contrast 
with these two contributions, the other two energy contributions due to the formation of K ion 
layers and to the deformation of FeSe layers are much smaller. Therefore, we conclude that 
the former two contributions dominate the energy change during the K intercalation. 

To understand the tendency for appearance of Fe vacancy in the K-intercalated iron 
selenides, we go further to consider the energy change AEre vacancy (= AE; - AE; - ufe, Where 


AE; = Exxre1.948e2 - Exorei.94se2 - XEx is the formation energy per unit cell in the Fe deficiency 


structure and ur, is the Fe chemical potential.) by removal of an Fe atom from the unit lattice 
of K,Fe2Se and its dependence on the K content. It is found that AER vacancy only fluctuates 
above and below zero when x < 0.6 as shown in Figure 2(a). When x > 0.6, AEre vacancy 
rapidly drops, revealing that Fe vacancies are favored at high levels of K intercalation. This is 
easily understood since K has a smaller electronegativity than Fe, when the doped electrons 
due to K is over a limit that Fe-Se bond can accommodate, the surplus electrons will prefer to 
transfer into Fe^*. In this way, Fe is repelled out of the lattice in the form of element. The final 
compound with Fe vacancy is more energetically favored state. The limit is around at x = 0.6. 
The K2Fe4Ses phase contains 20% Fe vacancy corresponding to x = 0.8 is a manifestation of 
this argument. This tendency of appearance of Fe vacancy with K content is supported by the 
experimental results Ea 

Besides the above considerations on the formation energy, the influence of lattice 
dynamics on the structural stability should also be accounted for. This is, in particular, true for 
the compounds with low levels of intercalated K, which can induce lattice instability. Figure 
2(b) shows the density of phonon states (DOPS) of Ko2Fe2Se2, where negative frequencies 
appear, meaning the structure is unstable. Partial DOPS further indicates that it is mainly 
because of the considerable amount of highly unstable Se atoms, which are unbonded due to 
the absence of K nearby. The phonon density, however, exhibits no imaginary mode for all 
other compounds with x > 0.25. As examples, Figure 2(c) shows the DOPS for Ko»sFe5Se» 
and Ko.sFe2Se2. The above results demonstrate that the phases of K,Fe2Se2 free of Fe vacancy 
are both energetically and dynamically favored in the K intercalation range of 0.25 < x < 0.6. 

Now we are focused on the trend of structural change upon K intercalation in the range 
of 0.25 < x < 0.6. Both electronic and size effects of K will be taken into accounts. In order to 
better understand the role of electron doping in the structural change, we charge KoFe»Se» 


with various electron concentrations, which allows us to explore the effect of electron doping 


while excludes the size influence of K. As shown in Figure 3(a) and (b), the changes in the 
lattice constant a, the Fe-Se bond and the Se-Fe-Se angle clearly indicate that FeSe layers are 
stretched along the ab plane upon electron doping alone. (Note the results of highly electron- 
doped are extrapolated from the trend of positively-charged ones simply because the stable 
negatively-charged KoFe2Se2 cannot always be obtained during the iterative calculations.) For 
a more realistic case, other factors such as the lattice mismatch and the Coulomb attraction 
must be included. Figure 3(c) shows the lattice parameter variations with K intercalation, 
which accounts for all these contributions together with the electron doping. We see that the 
overall effect of K intercalation, increases the lattice constant a, stretches FeSe layers and 
reduces the lattice constant c. For instance, the lattice constant a of Ko»sFe»Se» and Ko sFe>Se> 
expands from 3.65 Å to 3.69 Å and c shrinks from 14.42 Å to 13.92 Å. The predicted lattice 
constant c of Ko 25Fe2Se2, 14.42 Å, agrees with 14.28(4) Å of Ko3FezSe2 having T. of 44 K BH 
And the predicted lattice constants a and c for Ko sFe2Se>, 3.69 À and 13.92 Å, are consistent 
with the periods observed using scanning tunneling microscope along [110] and [001], 5.5 Å 
(V2a) and 14.1 Å P9, respectively considering the calculation accuracy. It is worth noting that 
the similar trend of the changes in lattice parameters was observed in other electron-doped 
ThCr5Sr structures, such as KFe>As> 9:44 and CaFe;As; Å. 

Furthermore, we explore the temperature-dependent stability of structures in the K 
region of interest. To this end, we perform the molecular dynamics (MD) simulations on 
Ko.25Fe2Se2 and Ko.5FeaSe2. Both structures can survive for at least 1 picosecond using 2V2 x 
2V2 x 1 supercells at temperatures up to 500 K, indicating that they are stable at this 
temperature 46] It should be worthy to note that the atom displacements of Ko.25Fe2Se2 are 
larger than those of Ko.sFe2Se2, suggesting that Ko.2sFe2Se2 is less stable. The variations of 
free energy with temperature for the two phases are also calculated and shown in Figure 3(d). 


The free energy is more favorable for Ko sFe>Se> at temperatures above 268 K, consistent with 


the result obtained by MD. Despite their variations in stability, both of them have the 
electronic structures similar to KFe2Se2 (see Figure S3), which suggests that they should have 
the similar properties to KFe2Se2. Considering stoichiometry, formation energy, stability and 
electronic structures, Ko.2sFe2Se2 and KosFe2Se2 (see Figure S2 for their structures) are 
proposed to responsible for the observed superconductivity at 44 K and at 30 K BN 
respectively. The relative stability difference could also explain the difficulty to obtain the 44 
K phase. 

Based on the results presented above, we schematize a phase diagram for the K,Fe2.ySe2 
system in Figure 4. In the K-rich portion with x > 0.6, phases with Fe vacancies tend to exist, 
agreeing well with the observed antiferromagnetic KoFe4Ses phase. As for the low level K- 
intercalated compounds (x < 0.25), there has been no report on the synthesis of free-standing 
K,Fe2Se2. We note that the two identified SC phases Ko3Fe2Se2 and KoFe2Se2 lie in the 
region of 0.25 < x < 0.6 P". T, is 44 K for the former phase and 30 K for the later one, 
suggesting that T. is dependent on K content or the doped electron concentrations. Although 
experimentally observed K contents in SC phases by far are discrete, recent reports about 


carrier concentration tuning of T. from 30 to above 40 K in FeSe thin flakes [47, 48] 


implies that 
their variation of T. with carrier concentration can be continuous, similar to the cases in other 
high-temperature superconductors. Further improvement of T. can be expected considering 


[49-54] Therefore, we infer that 


the experimental progresses achieved in single-layer FeSe films 
the 30 K phase is electron overdoped and the 44 K phase also might not be optimally tuned. 
Moreover, in alkali-metal-intercalated FeSe compounds prepared at low temperatures, the 
synergic effect of NH3, NH», or C2H4(NH2) along with alkali metal can stabilize the 
structures P^ 5^ 55! These offer us an effective strategy to raise T. of bulk iron-selenide-based 


superconductors by controlling carrier doping while stabilizing the structures by intercalations 


or effect of substrates. At the moment, tremendous efforts are needed toward this goal. 


Because of the calculation limit and accuracy, we do not consider the slight off-stoichiometric 
cases, say, the Fe vacancy concentration is less than 3 at.%. It should be pointed out that such 
slight off-stoichiometry is tolerated in K;Fe»Se» just like many other materials, which may be 
the reason that superconductivity was observed in previous reports of refs. [33, 40]. 

In conclusion, we carefully investigate the energetic change and structural evolution of 
K,Fe>.,Se> as a function of intercalated K content using first-principles calculations. Two 
factors dominating the formation of the phases and the structural evolution are confirmed. 
One is due to the accumulation of negative charge in FeSe layers; the other is due to Coulomb 
attraction between K ion layers and negatively charged FeSe layers. The structural evolution 
of this series of phases is summarized: at 0.25 < x < 0.6, the bct lattice is stable and Fe is 
favored to full occupy its site; at x > 0.6, FeSe layers tend to exclude Fe atoms and create Fe 
vacancies; and at x < 0.25, the structure will collapse for the dynamic instability. A schematic 
phase diagram is constructed accordingly and the possible route to further improve T. is 
suggested. The phases responding to the observed superconductivity are proposed to be 
Ko25Fe2Se2 and KosFe2Se2 in terms of stoichiometry, formation energy, stability and 
electronic structures. Though based on the study of K,Fe>.,Se>, our results should be 
meaningful to understand the SC and related phases in metal-intercalated iron selenides and 


other similar SC systems. 
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Figure 1. (a) The schematic of the K intercalation process. It does not represent a certain 
structure. (b) The total formation energy and the energy change due to the formation of K ion 
layers, the electron doping in FeSe layers, the deformation of FeSe layers and the Coulomb 


attraction between K ion layers and FeSe layers as a function of K content, respectively. 
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Figure 2. (a) The relative variation of formation energy change between K,Fe;s4Se; and 
K,Fe2Se, as a function of K content. As the absolute value will be determined by Fe chemical 
potential (AEFe vacancy = AE; - AE, - Ure), the change of formation energy between 
Ko.125Fe;.94Se? and Ko.125Fe2Se2 is set to zero to show the relative values for comparison. (b) 


Total and partial DOPS of Ko?FeoSe». (c) DOPS of Ko.25Fe2Se2 and Ko.sFe2Sez. 
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Figure 3. The lattice constant a (a), the Fe-Se bond length and the Se-Fe-Se bond angle (see 
details in the inset) along FeSe layers (b) as a function of charge per unit cell. Positive charge 
means hole-doped and negative one being electron-doped. Dashed lines represent the results 
of linear fitting. (c) The variation of the lattice constants a (pink rings) and c (orange stars) 
with the K content in K,Fe2Se2. (d) The free energy of Ko5sFe;Se; and Ko.sFe2Se2 as a 


function of temperature. 
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1. Simplification of K,Fe;Se; superstructures 


The example about the minor energy difference of superstructures due to different 
arrangements of K ions along c axis is shown in Figure S1(a) and S1(b). The example 
about supercells with different lattice constants in ab plane have higher AE; compared 


to those with the same lattice constants is shown in Figure S1(c) and S1(d). 


2. The proposed structures of K,F&Se; (x = 0.3 and 0.6) 


The possible structures of K,Fe2Se2 with uniform K distribution are proposed, 
Ko25sFe2Se2 and KosFe2Se2. For x = 0.3, Ko25sFe2Se2 with K ion layers in a 
superstructure of 2 x 2 is proposed and the arrangement of K ions is shown in Figure 


S2(a). In Figure S2(b), K,Fe2Se2 (x = 0.6) is with the superstructure V2 x N2 of K ion 


layers, which is actually with composition of Ko,sFe2Se2. This may be the structure 
mentioned in previous reports 121 AE; of Ko5sFe5Se» (-0.87 eV/unit cell) is higher 
than that of Ko sFe2Se2 (-1.62 eV/unit cell), implying that Ko.25Fe2Se2 should be less 


stable than Ko 5Fe>Se> at the ground state. 


3. The electronic structures of K,Fe2Se, (x ~ 0.3 and 0.6) 


Calculations of total and partial density of states (DOS) are performed to 
investigate the effects of superstructures on the properties. In Figure S3(a) and S3(b), 
Fe 3d and Se 4p orbital hybridization is located between -7 to -3 eV, which is 
responsible for the formation of FeSe layers in both Ko.sFe>Se> and Ko.sFe2Se2. Part 
of Fe 3d orbital is located around the Fermi level and decides the properties of both 
structures. As shown in Figure S3(c), the main difference among total DOS of 
Ko.25Fe2Se2, KosFe2Se2, and KFe2Se2 is the locations of their Fermi levels. If the 
Fermi level of KFe>Se» is set to zero, the values of Fermi levels for Ky>5Fe>Se> and 
Ko.sFe2Se2 are about -260 and -200 meV, respectively. So the properties of Ko.25Fe2Se2 
and Ko.5Fe2Se2 should be basically the same with that of KFe2Se2, while their carrier 
concentrations are different so that they may have different transition temperatures in 


terms of superconductivity. 


4. Methods 


All the first-principles calculations employ generalized gradient approximation 


based on density functional theory in the form of Perdew-Burke-Ernzerhof function x 


for exchange-correlation potential, which are performed in the Cambridge Serial Total 
Energy Package |“. The self-consistent field method is used with a tolerance of 5.0 x 
10” eV/atom, which is in conjunction with plane-wave basis sets of cutoff energy of 


330 eV and ultrasoft pseudopotentials in reciprocal space Dl 


According to 
Monkhorst-Pack special k-point scheme 6! the first Brillouin zone is sampled with 
grid spacing of 0.04 ÀT. Fully optimization of the atomic positions and lattice 
parameters of all compounds is done until the remanent Hellmann-Feynman forces on 
all components are less than 0.01 eV/Å. The finite displacement method is used to 
calculate density of phonon states. In order to examine the dynamic stability of the 
structures, we perform first-principles finite temperature molecular dynamics 
simulations with time steps of 10 femtosecond under pressure of 0.1 MPa. Free 
energies are calculated as a function of temperature using the quasiharmonic 
approximation based on phonon dispersion ". 

The initial model of KFe>Se» is isostructural to BaFe>As» [8], Superstructures due 
to different arrangements of K ions along c axis are not considered for the minor 
energy difference. For example, the difference of AE; between two kinds of 1x 1 x 2 
supercells with composition of Ko,sFe2Se2 [see Figure S1(a) and S1(b)], is only 0.4 
meV/unit cell. Superstructures containing one K layer and one FeSe layer due to 
different arrangements of K ions perpendicular to c axis are the main concerns. The 
AE; of 4 x I x 1 supercell with composition of Ko,sFe2Se2 [see Figure S1(c)] is 26 


meV/unit cell higher than that of 2 x 2 x 1 supercell with the same composition [see 


Figure S1(d)]. Hence, supercells with different lattice constants in ab plane, like 3 x I 


x 1,3 x 2 x I and so on, are not discussed in this work for their higher AE; compared 
to those with the same lattice constants in ab plane. Supercells V2 x V2 x 1,2x2x1, 
V5 x N5 x I and 2V2 x 2V2 x 1 are built accordingly with uniform K distribution as 
the Coulomb repulsion between K ions dominates K distribution. To study the 
formation of Fe vacancy, a 2V2 x 2V2 x 1 supercell (a = b > 10.3 A and c > 13.5 Å), 
which is large enough to accommodate a single Fe vacancy while keeping the 
tetragonal symmetry, is therefore chosen. In this case, the Fe content corresponds to y 
— 0.06. The above structures are only designed to obtain different K contents in 
K,Fe».,Se» (y = 0 or 0.06), which shall not represent the real ones. As K ion ordering 
cannot guarantee the tetragonal symmetry in the structures, the obtained lattice 
constants or bond lengths will be the average values. The bond angles are obtained 
using the average coordinates of Fe and Se atoms. Besides, when we calculate AF, 
the total energy of KoFe2Se2 is used instead of that of FeSe so as to offset the 


additional energy, tens of meV per unit cell induced by the change of cell volume. 
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Figure S1 (a) (b) The structures of two kinds of 1 x 1 x 2 supercells with composition 
of Ko,sFe2Se2. The former structure always has a vacancy between two K ions while 
the K ions in the latter appear in pairs along c axis. (c) The structure of 4 x 1 x I 
supercell with composition of Ko.5Fe2Se2. (d) The structure of 2 x 2 x 1 supercell with 
composition of Kos5Fe2Se2. Violet (largest), grey (smallest) and orange spheres 


represent K, Fe and Se atoms, respectively. 


LI 
o 
e 
o 
1 e 
> © 
e 
o. ©. o- o 
‘ec ED 
©. ©. © `o 
w e o e o 
e - e - e - e 
We 


Figure S2 The possible K ion arrangements in Kọo.25Fe2Se2 (a) and Ko,sFe2Se2 (b). 


Violet big spheres represent K ions and orange small spheres represent vacancies. 
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Figure S3 (a) (b) The partial DOS of Fe 3d and Se 4p in Ko25Fe2Se2 and Ko. 5Fe2Se2. 
The Fermi level is set to zero. (c) The comparison of total DOS details around the 
Fermi levels of Ko.25Fe2Se2, Ko, sFe>Se, and KFe2Se2. The Fermi level of KFe,Se; is 


set to zero and others are set accordingly. 


